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Abstract. We present ID radiative transfer modelling of the envelopes of a sample of 18 low-mass protostars and 
pre-stellar cores with the aim of setting up realistic physical models, for use in a chemical description of the 
sources. The density and temperature profiles of the envelopes are constrained from their radial profiles obtained 
from SCUBA maps at 450 and 850 (jm and from measurements of the source fluxes ranging from 60 /im to 1.3 
mm. The densities of the envelopes within ~10000 AU can be described by single power-laws p oc r~ a for the 
class and I sources with a ranging from 1.3 to 1.9, with typical uncertainties of ±0.2. Four sources have flatter 
profiles, either due to asymmetries or to the presence of an outer constant density region. No significant difference 
is found between class and I sources. The power-law fits fail for the pre-stellar cores, supporting recent results 
that such cores do not have a central source of heating. The derived physical models are used as input for Monte 
Carlo modelling of submillimeter C ls O and C 17 emission. It is found that class I objects typically show CO 
abundances close to those found in local molecular clouds, but that class sources and pre-stellar cores show 
lower abundances by almost an order of magnitude implying that significant depletion occurs for the early phases 
of star formation. While the 2-1 and 3-2 isotopic lines can be fitted using a constant fractional CO abundance 
throughout the envelope, the 1-0 lines are significantly underestimated, possibly due to contribution of ambient 
molecular cloud material to the observed emission. The difference between the class and I objects may be related 
to the properties of the CO ices. 
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1. Introduction 

In the earliest, deeply-embedded stage a low- mass proto- 
star is surrounded by a collapsing envelope and a circum- 
stellar disk through which material is accreted onto the 
central star, while the envelope is dissipated simultane- 
ously through the action of the powerful jets and outflows 
driven by the young star. Traditionally, young stellar ob- 
jects (YSOs) have been classified according to their spec- 
tral energy distributions (SEDs) in the class I-III scheme 
(Lada, 1987; Adams et al., 1987) describing the evolu- 



of |Myers et al. (1983| ) and |Bcnson fc Myers (1989| ) are 



good candidates for pre-stellar cores, i.e., dense gas cores 
that are on the brink of collapse and so leading to the 
class and I phases. 



tion of YSOs from the young class I sources to the more 
evolved pre-main sequence class III sources. This classifi- 



The |Shu (1977| ) model predicts that the outer parts 
of the envelope follow a p oc r~ 2 density profile similar 
to the solution of an isothermal sphere (Larson, 1969), 
while within the collapse radius, which is determined by 
the sound speed in the envelope material and the time 
since the onset of the collapse, the density tends to flatten 
nearing a p oc r profile in the innermost parts. This 
model has subsequently been refined to include for exam- 



cation scheme was further expanded by Andre et al. (1993) 
to include sources that mainly radiate at submillimeter 
wavelengths (i.e., with high ratios of their submillimeter 
and bolometric luminosities, -Zj S ubmm/£boi) and it was sug- 
gested that these so-called class sources correspond to 
the youngest deeply embedded protostars. Even earlier in 
this picture of low-mass star formation, the starless cores 
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pie rotational flattening (Terebey et al., 1984) where such 
effects are described as a perturbation to the Shu (1977) 
solution. 

An open question about the properties of YSOs in the 
earliest protostellar stage is how the structure of the en- 
velope reflects the initial protostellar collapse and how it 
will affect the subsequent evolution of the protostar, for 
example in defining the properties of the circumstellar disk 
from which planets may be formed later. One of the pos- 
sible shortcomings of the Shu-model is the adopted, con- 
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stant, accretion rate. Foster & Chevalier (1993) performed 
hydrodynamical simulations of the stages before the pro- 
tostellar collapse and found that the structure of the core 
at the collapse initiating point is highly dependent on the 
initial conditions; only in the case where a large ratio ex- 
ists between the radius of an outer envelope with a flat 
density profile and an inner core with a steep density pro- 
file, will the core evolve to reproduce the conditions in the 



Shu model. In an analytical study, Henriksen et al. (1997) 
suggested that the accretion history of protostars could 
be divided into two phases for cores with a flat inner den- 
sity profile: a violent early phase with high accretion rates 
(corresponding to the class phase) that declines until 
a phase with mass accretion rates similar to the predic- 
tions in the Shu- model is reached (class I objects), i.e., a 
distinction between class and I objects based on ages. 
Whether t his is indeed the case has recently been ques- 
tioned by Jayawardhana et al. (2001 ), who instead sug- 
gest that both class and I objects are protostellar in 
nature, but just associated with environments of differ- 
ent physical properties, with the class objects in more 
dense environments leading to the higher accretion rates 
observed towards these sources. 

The chemical composition of the envelope may be 
an alternative tracer of the evolution. Indeed, for high- 
mass YSOs, combined infrared and submillimcter data 
have shown systematic heating trends reflected in the ice 
spectra, gas/ice ratios and gas-phase abundances (e.g 



Gerakines et al., 1999; Boogert et al., 2000; van der Tak 



et al., !000a). One of the prime motivations for this work 



is to extend similar chemical studies to low-mass objects, 



and extensive (sub-)millimeter line data for a sample of 
such sources are being collected at various telescopes, 



ing density distribution is also a power-law (i.e., of the 
type p ex r~ a ), one can then derive a relationship between 
the radial intensity observed in continuum images and the 
envelope radius, which will also be a power-law with an 
exponent depending on the power-law indices of the den- 
sity and temperature distributions. Both groups find that 
the data sets are consistent with a in the range 1.5-2.5 
in agreement with previous results and the model predic- 
tions. However, as both groups also notice, in the case 
where the assumption about an optically thin envelope 
breaks down, the temperature distribution and so the de- 
rived density distribution may not be correctly described 
in this approach. 

To further explore these properties of the protostellar 
envelopes we have undertaken full ID radiative transfer 
modelling of a sample of protostars and pre-stellar cores 
( see Sect. using th e radiative transfer code DUSTY 
(Ivezic & Elitzur, 1997). Assuming power-law density dis- 
tributions we solve for the temperature distribution and 
constrain the physical parameters of the envelopes by com- 
parison of the results from the modelling to SCUBA im- 
ages of the individual sources and their spectral energy 
distributions (SEDs) using a rigorous \ 2 method. Besides 
giving a description of the physical properties of low-mass 
protostellar envelopes, the derived density and temper- 
ature profiles are essential as input for detailed chem- 
ical modelling of molecules observed towards these ob- 
jects. Also, a good description of the envelope structure is 
needed to constrain the properties of the disks in the em- 
bedded phase (e. g., Keene fc Masson, 199C ; Hogcrheijde 
et al., 1998|Jl999| ; [Looney et al., 2000j 



In Sect. Hour sample of sources is presented and the re- 



which can be complemented by future SIRTF infrared 
data. 

In order to address these issues the physical parameters 
within the envelope, in particular the density and temper- 
ature profiles and the velocity field, are needed. The first 
two can be obtained through modelling of the dust con- 
tinuum emission observed towards the sources, while ob- 
servations of molecules like CO and CS can trace the gas 
component and velocity field. At the densities observed in 



duction and calibration briefly discussed (see also Schoici 
et al., 2002 ). In Sect. || the modelling of the sources is 



described and the derived envelope parameters presented. 
The prop erties of the individual sources are described in 
Sect. 3.4. In Sect. the implication of these results are 



discussed and compared to other work done in this field. 
The results of the continuum modelling will be used in a 
later paper as physical input for detailed radiative transfer 
modelling of molecular line emission for the class objects 



the inner parts of the envelopes of YSOs it is reasonable to 



expect gas and dust coupling, which is usually expressed 
by the canonical dust-to-gas ratio of 1:100 and the as- 
sumption that the dust and gas temperatures are similar. 
Therefore a physical model for the envelopes derived on 
the basis of the dust emission can be used as input for 
modelling of the abundances of the various molecules. 

Recently, |Shirley et al. (2000| ) and |Motte fc Andre 
(200l|^have undertaken surveys of the continuum emis- 
sion of low-mass protostars - using respectively SCUBA 
(at 450 and 850 fim) and the IRAM bolometer at 1.3 
mm. Both groups analyzed the radial intensity profiles 
(or brightness profiles) for the individual sources, assum- 
ing that the envelopes are optically thin, in which case the 
temperature follows a power-law dependence with radius 
in the Rayleigh- Jeans limit. Assuming that the under ly- 



as has been done for class I objects ( Hogerheijde et al 
1998| ) and high-mass YSOs ( |van der Tak et al., 2000b| ) and 



as presented for the low-mass class object IRAS 16293- 
2422 flCeccarelli et al., 20004 |b] |Schoicr et al., 2002| ). In 
Sect. |5j, the first results of this radiative transfer analysis 
for C O and C 17 is presented. 



2. Data, reduction and calibration 

2.1. The sample 

The class sources in the sample have been chosen from 



the list of Andre et al. (2000), with the requirement that 
they should have a distance of less than 450 pc, luminos- 
ity less than 50 L Q and be visible from the JCMT. In 



addition, CB244 from the Shirley et al. (200C) sample was 
included. These objects were supplemented with two pre- 
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stellar cores L1689B and LI 544, also from Shirley et al 
(20001 )] and two class I objects, L14 89 and TMR1 taken 
from Hogerheijde et al. (1997, 1998). To enlarge the class 
I sample, L1551-I5, TMC1A and TMC1 were included as 
well. The physical properties for these three sources were 
modelled using the same approach as the remainder of the 
sources, based on SCUBA archive data. They were, how- 
ever, not included in the JCMT line survey, so the line 
modelling (Sect. was mainly based on data presented 
in the literature, in particular Hogerheijde et al. (1998) 



were found between the individual maps. Two sources, 
L1157 and CB244 only had usable data at 850 /jm (see 
also Shirley et al., 200"o| ) , so supplementary data for these 



and Ladd et al. (1998) 



For the class objects we have adopted luminosities 
and distances from |Andrc et al. (200C|), for the class I 



objects the values from Motte & Andre (2001) and for 
the prc-stellar cores and CB244 distances and luminosities 
from Shirley et al. (2000). There are a few exceptions, 



however: for the objects related to the Perseus region we 
assume a distance of 220 pc and scale the luminosities 
from Andre et al. (200"o| ) accordingly, while a distance of 
325 pc is assumed for L1157 as in IShirlcy et al. (2000| ). 
The sample is summarized in Table 1. The class object 
IRAS 16293-2422 treated in [Schoicr et al. (2002[ ) has been 



included for comparison here as well. 
2.2. Submillimeter continuum data 

Archive data obtained from the Submillimctre Common- 
User Bolometer Array, (SCUBA), on the James Clerk 
Maxwell TelescopeQ (JCMT), on Mauna Kea, Hawaii were 
adopted as the basis for the analysis. Using the 64 bolome- 
ter array in jiggle mode, it is possible to map a hexagonal 
region with a size of approximately 2.3' simultaneously at, 
e.g., 450 [im and 850 fim. It is also possible to combine 
jiggle maps with various offsets to cover a larger region. 
To perform the initial reduction of the data, the package 



SURF (Icnncss & Lightfoot, 1997) was used following the 
description in Sandcll (1997|) . The individual maps were 
extinction corrected with measurements of the sky opac- 
ity r obtained at the Caltech Sub millimeter Observatory 
(CSO) and using the relations from Archibald et al. (2000|) 
to convert the CSO 225 GHz opacity to estimates for the 
sky opacity at 450 /Ltm and 850 fim. The sky opacities 
can also be estimated using sky dips, and in cases where 
these were obtained, the two methods agreed well. Most of 
the sources were observed in the course of more than one 
program and on multiple days, so wherever possible avail- 
able data obtained close in time were used, coadding the 
images to maximize the signal-to-noise and field covered. 
In the coadding, it is possible to correct for variations in 
the pointing by introducing a shift for each image found 
by, e.g., fitting gaussians to the central source. We chose, 
however, not to do this, because only minor corrections 



1 The JCMT is operated by the Joint Astronomy Centre 
in Hilo, Hawaii on behalf of the parent organisations: the 
Particle Physics and Astronomy Research Council in the 
United Kingdom, the National Research Council of Canada 
and the Netherlands Organization for Scientific Research. 



two sources were obtained in October 2001 at 450 /im (see 
Sect. fD|). 

For each source the flux scale was calibrated using 
available data for one of the standard calibrators, either 
a planet or a strong submillimeter source like CRL618. 
From the calibrated maps the total integrated fluxes were 
derived and the ID brightness profiles were extracted by 
measuring the flux in annuli around the peak flux. The 
annuli were chosen with radii of half the beam (4.5" for 
the 450 ^m data and 7.5" for the 850 ^m data) so that 
a reasonable noise-level is obtained, while still making 
the annuli narrow enough to get information about the 
source structure without oversampling the data. Actually 
the spread in the fluxes measured for the points in each 
annulus due to instrumental and calibration noise was neg- 
ligible compared to the spread due to (1) the gradient in 
brightness across each annuli, and (2) deviations from cir- 
cular structure of the sources. One problem in extract- 
ing the brightness profiles was presented by cases where 
nearby companions were contributing significantly when 
complete circular annuli were constructed - the most ex- 
treme example being N1333-I4 with two close protostars. 
In these cases emission from "secondary" components was 
blocked out by simply not including data-points in the di- 
rection of these closeby sources when calculating the mean 
flux in each annulus. 

2.3. SCUBA observations of 11157 and CB244 



The observations of LI 157 and CB244 were obtained on 
October 9th, 2001. Calibrations were performed by ob- 
serving Mars and the secondary calibrator, CRL2688, 
immediately before the observations. Skydips were ob- 
tained immediately before the series of observations (all 
obtained within 3 hours) giving values for the sky opac- 
ity of T450 = 1.2 and t^q = 0.23, which agree well with 
the sky opacities estimated at the CSO during that night. 
From gaussian fits to the central source the conversion 
factor from the V onto the Jy beam -1 scale (C\) was es- 
timated and is summarized in Table ^| together with the 
beam size (9 m b also estimated from the gaussian fit to the 
calibration source. For Mars the estimate of the beam size 
was obtained by deconvolution with the finite extent of the 
planet, while C RL2688 was assumed to be a point source 
( ISandcll, 19941 ). 

The derived parameters for L1157 and CB244 are given 
in Table ||. Images of the two sources at the two SCUBA 
wavelengths are presented in Fig. |l|. As seen from the fig- 
ure, both sources are quite circular with only a small de- 
gree of extended emission. C omparison with the 850 /im 
data of [Shirley et al. (2000 ) for the 40" aperture shows 
that the fluxes agree well within the 20% uncertainty as- 
sumed for the calibration. 
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Table 1. Sample of sources. 





a (2000) 


5(2000) 


Tboi 




d 


Other names 


Type 




(hh mm ss) 


(dd mm ss) 


(K) 




(pc) 






L1448-I2 


03 25 22.4 


+ 30 45 12 


60 


3 


220 


03222+3034 


Class 


L1448-C 


03 25 38.8 


+ 30 44 05 


54 


5 


220 


L1448-MM 




N 1333-12 


03 28 55.4 


+ 31 14 35 


50 


16 


220 


03258+3104/SVS19 




N 1333-14 A 


03 29 10.3 


+ 31 13 31 


34 


6 


220 






N1333-I4B 


03 29 12.0 


+ 31 13 09 


36 


6 


220 






L1527 


04 39 53.9 


+ 26 03 10 


36 


2 


140 


04368+2557 




VLA1623 


16 26 26.4 


- 24 24 30 


35 


1 


160 






L483 


18 17 29.8 


- 04 39 38 


52 


9 


200 


18148-0440 




L723 


19 17 53.7 


+ 19 12 20 


47 


3 


300 


19156+1906 




L1157 


20 39 06.2 


+ 68 02 22 


42 


6 


325 


20386+6751 




CB244 


23 25 46.7 


+ 74 17 37 


56 


1 


180 


23238+7401/L1262 




IRAS 16293-2422 


16 32 22.7 


- 24 28 32 


43 


27 


160 




(a) 


L1489 


04 04 43.0 


+ 26 18 57 


238 


3.7 


140 


04016+2610 


Class I 


TMR1 


04 39 13.7 


+ 25 53 21 


144 


3.7 


140 


04361+2547 




L1551-I5 


04 31 34.1 


+ 18 08 05 


97 


28 


140 


04287+1801 


W 


TMC1A 


04 39 34.9 


+ 25 41 45 


172 


2.2 


140 


04365+2535 


W 


TMC1 


04 41 12.4 


+ 25 46 36 


139 


0.66 


140 


04381+2540 


(b) 


LI 544 


05 04 17.2 


+ 25 10 44 


18 


1 


140 




Pre-stellar 


L1689B 


16 34 49.1 


- 24 37 55 


18 


0.2 


160 







Notes: (a) Class object treated in |Schoicr ct al. [2002}. (h) Cl ass I object not included in the JCMT line survey but with CO 
observations from |Hogcrhcijdc et aL*7l998^ Laddet al. (1998j ). 



Table 2. Summary of the calibration for the October 2001 
data. 





A 










Mars 


450 


/Mm 


1.2 


395.1 


8.9" 


CRL2688 


450 


/im 


1.2 


310.5 


9.1" 


Mars 


850 


fim 


0.23 


297.1 


15.3" 


CRL2688 


850 




0.23 


259.6 


15.3" 



Notes: ' a 'Sky opacity. ' 'Conversion factor from V to Jy 
beam -1 scale. ' c 'Beam size (HPBW). 

Table 3. Results for the CB244 and LI 157 submillimeter 
emission. 



CB244 L1157 
450 fim 850 fim 450 ^m 850 /jm 



■ftp cak ^ 


3.14 


0.591 


8.62 


1.72 


p 0) 

-t+40" 


14.2 


1.28 


22.2 


2.74 


F Ijl2 o» C6) 


61.4 


3.11 


69.2 


5.87 


TP . ( a ) 
* noise 


0.43 


0.087 


0.29 


0.074 



Notes: (a) Peak flux and RMS noise in Jy beam . 1 'Integrated 
flux in 40" and 120" apertures respectively in Jy. 

2.4. Line data 

CO line data were obtained with the JCMT in May 
and August 2001, the fnstituto de Radio Astronomica 
Milimetrica (IRAM) 30m telescope in November 2001 and 
the Onsala Space Observatory 20m telescope in March 
2002, complementing data from the JCMT archive. Our 
own observations were performed in beam switching mode 



using a switch of 180" in declination - except for the 
sources in NGC1333, where position switching towards an 
emission free reference position was used. A more detailed 
description of the JCMT and the heterodyne receivers can 
be found on the JCMT homepageP| Where archive data 
were available for one line from several different projects, 
the data belonging to each observing program were re- 
duced individually and the results compared giving an es- 
timate of the calibration uncertainty of the data of 20%. 
The integrated line intensities were found by fitting gaus- 
sians to the main line. In some cases outflow or secondary 
components were apparent in the line profiles leading to 
two gaussian fits. For the C 17 J = 1 — and J = 2 — 1 
lines, the hyperfine splitting were apparent, giving rise to 
two separate lines for the J = 1 — transition separated by 
about 5 km s , while the J — 2 — 1 main hyperfine lines 
are split by less (0.5 km s _1 ) giving rise in some cases to 
line asymmetries. In these cases the quoted line intensities 
are the total intensity including all hyperfine lines. 

The integrated line data were brought from the an- 
tenna temperature scale T£ to the main-beam brightness 
scale T m b by dividing by the main-beam brightness effi- 
ciency ?7 m b taken to be 0.69 for data obtained using the 
JCMT A band receivers (210-270 GHz; the J = 2-1 tran- 
sitions) and 0.59-0.63 for respectively the old B3i (before 
December 1996) and new B3 receivers (330-370 GHz; the 
,7 = 3 — 2 transitions). For the IRAM 30m observations 
beam efficiencies B c g of 0.74 and 0.54 and forward efficien- 
cies F e s of 0.95 and 0.91 were adopted for respectively the 
C 17 J = 1 — and J = 2 — 1 lines, which corresponds 

2 tittp : //www . j ach . hawaii . edu/ JACpublic/ JCMT/ 
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CB244 
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850 fim 
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O 
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RA offset ["] RA offset ["] 

Fig. 1. SCUBA images of CB244 and L1157 at 450 and 850 /im. The contours indicate the intensity corresponding to 
2cr, 4cr, etc. with a being the RMS noise given for each source and wavelength in Table |[ 



to main-beam brightness efficiencies (?? m b = B c ff/F e g) of 
respectively 0.78 and 0.59. For the Onsala 20m telescope 



ulated dust grains with thin ice mantles at a density of 
10 6 cm -3 were adopted. These were found by van 



«H 2 



r/ m b = 3.43 was adopted for the C 18 and G r '0 J = 1 — dcr Tak ct al. (1999] ) to be the only dust opacities that 



lines. The relevant beam sizes for the JCMT are 21" and 
14" at respectively 220 and 330 GHz, for the IRAM 30m, 
22" and 11" at respectively 112 and 224 GHz and for 
the Onsala 20m, 33". The velocity resolution ranged from 
0.1-0.3 km s" 1 for the JCMT data and were 0.05 and 
0.1 km s _1 for the observations of respectively the C 17 
J = 1 - and J = 2 - 1 transitions at the IRAM 30m. 
The line properties are summarized later in Sect. ^[ 

3. Continuum modelling 

3.1. Input 

To model the physical properties of the envelopes around 
these sources the ID radiative transfer code DUSTY 
(Ivezic et al., 1999 ) was used.PI T he dust grain opacities 



from Ossenkopf & Hcnning (1994) corresponding to coag 



DUSTY is publically available from the homepage at 
http : //www.pa.uky . edu/~moshe/dusty/ 



could reproduce the "standard" dust-to-gas mass ratio 
of 1:100 by comparison to C 17 measurements for warm 
high-mass YSOs where CO is not depleted. 

Using a power-law to describe the density leaves five 
parameters to fit as summarized in Table |[ Not all five 
parameters are independent, however: the temperature at 
the inner boundary, Xi, determines the inner radius of 
the envelope, r\, through the luminosity of the source. 
If the outer radius of the envelope ri is expected to be 
constant, Y = r-ilr\ will depend on the value of r±, i.e., 
T\. The results are, however, not expected to depend on 
t i if it is chosen small enough, since the beam size does 
not resolve the inner parts anyway. Therefore Tj is simply 
set to 250 K, a reasonable temperature considering the 
chemistry observed towards these sources. Also the tem- 
perature of the central star has to be fixed: a temperature 
of 5000 K is chosen. This temperature is of course mostly 
unknown for the embedded sources, but due to the optical 
thickness of the envelope most of the radiation from the 
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Table 4. Parameters for the DUSTY ID radiative transfer 
modelling of the protostellar envelopes. 

Param. Description 
Modelled: 

Y Ratio of the outer (7-2) to inner (n) radius 

Tloo Optical depth at 100 fim 

a Density power-law exponent 



Fixed: 
Ti 
T* 

Literature: 

d Distance 

Z/boi Luminosity 



Temperature at inner boundary (250 K) 
Temperature of star (5000 K) 



central star is anyway reprocessed by the dust and thus 
the temperature of the central star does not play a critical 
role, e.g., in the resulting SED. 

Although these parameters may not seem the most 
straightforward choice, one of the advantages of DUSTY 
is the scale-free nature allowing the user to run a large 
sample of models and then compare a number of YSOs to 
these models just by scaling with distance and luminosity 



as discussed in Ivczic et al. (1999) 



3.2. Output 

DUSTY provides fluxes at various wavelengths and bright- 
ness profiles for the sources, which are compared to the 
SCUBA data and flux measurements. Given the grid of 
models the best fit model can then be determined by 
calculating the x 2 -statistics for the SED and brightness 
profiles at 450 and 850 fim (xIedj X450 an d xiso respec- 
tively). In order to fully simulate the observations, the 
modelled brightness profiles are convolved with the exact 
beam as obtained from planet observations. Strictly speak- 
ing, the outer parts of the brightness profile also depend 
on the chopping of the telescope. The chopping along one 
axis does by nature not obey the spherical symmetry, so 
simulation of the chopping and comparing this with one 
dimensional modelling will not reflect the observations. 
Therefore in this analysis only the inner 60" of the bright- 
ness profiles arc considered, which (1) should be less sensi- 
tive to the typical 120" chop and (2) is typically above the 
background emission. For the flux measurements a relative 
uncertainty of 20% was assumed irrespective of what was 
given in the original reference, since some authors tend to 
give only statistical errors and do not include calibration 
or systematic errors. By assigning a relative uncertainty 
of 20% to all measurements each point is weighted equal 
but more weight is given to a given part of the SED if 
several independent measurements exist around a certain 
wavelength. Contour plots of the derived x 2 values for 
L483-mm are presented in Fig. ^|, while the actual fits to 
the brightness profiles and the SED for this source are 
shown in Fig. |^ and ^. 



In determining the best fit model each of the calcu- 
lated x 2 values are considered individually. The total % 2 
obtained by adding the Xsed: xlso an d X450 does not make 
sense in a strictly statistical way, since the observations 
going into these cases are not 100% independent. Another 
reason for not combining the values of Xsedi X450 an d 
Xs,50 mT, ° on e total x 2 is that the parameters constrained 
by the SED and brightness profiles are different. For ex- 
ample the brightness profiles provide good constraints on 
a as seen in a 2D contour (Y, a) plot of, e.g., X450 m Fig. |[ 
while these do not depend critically on the value of Tioo- 
The most characteristic feature of the x 2_va l ue s for the 
SEDs on the other hand is the band of possible models 
in contour plots for (a,rioo), giving an almost one-to-one 
correspondence for a best fit tioo for each value a. 

These features are actually easily understood: X450 an d 
x| 50 are the normalized profiles and should thus not de- 
pend directly on the value of Tioo- On the other hand since 
the peaks of the SEDs are typically found at wavelengths 
longer than 100 /zm, increasing tioo and thus the flux at 
this wavelength, require the best-fit SED to shift towards 
100 /im, i.e., with less material in the outer cool parts of 
the envelope, which can be obtained by a steeper value 
of a. The value of Y is less well constrained, mainly be- 
cause of its relation to the temperature at the inner radius 
(and through that the luminosity of the central source). 
As illustrated in Fig. |^, Y is constrained within a factor 
1.5-2.0 at the 2a level, but the question is how physical 
the outer boundary of the envelope is: is a sharp outer 
boundary expected or rather a soft transition as the den- 
sity and temperature in the envelope reaches that of the 
surrounding molecular cloud? In the first clear drop 

of the observed brightness profile should be seen compared 
with a model with a (sufficiently) large value of Y, e.g., 
corresponding to an outer temperature of 5-10 K, less 
than the temperature of a typical molecular cloud. In the 
other case, however, such a model will be able to trace 
the brightness profile all the way down to the noise limit. 
The modelling of the CO lines (see Sect. ||) indicates that 
significant ambient cloud material is present towards most 
sources, so the transition from the isolated protostar to the 
parental cloud is likely to be more complex than described 
by a single power-law. 

The features of the values of Xsedj X450 anc i Xsso P r0 " 
vide an obvious strategy for selecting the best fit mod- 
els: first the best fit value of a is selected on the basis 
of the brightness profiles and second the corresponding 
value of Tioo is selected from the Xsed contour plots. For 
a few sources there is not 100% overlap between the 450 
and 850 /im brightness profiles and it is not clear which 
brightness profile is better. The beam at 850 /zm is signif- 
icantly larger than that at 450 /xm (15" vs. 9") and even 
though the beam is taken into account explicitly, the sen- 
sitivity of the 850 fim. data to variations in the density 
profile must be lower as is seen from Fig. ||. On the other 
hand, the 450 /im data generally suffer from higher noise, 
so especially the weak emission from the envelope, which 
provides the better constraints on the outer parts of the 
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Fig. 2. x 2 contour plots for the modelling of L483-mm. In 
the four upper panels Tioo is fixed at 0.2, in the 4 middle 
panels a is fixed at 0.9 and while in the 4 lower panels 
Y is fixed at 1400. The solid (dark) contours indicate the 
confidence limits corresponding to la, 2a etc. 



envelope and thus the power-law exponent, will be more 
doubtful at this wavelength. The power-law slopes found 
from modelling the two brightness profiles agree, however, 
within the uncertainties (a ~ ±0.2). 
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Fig. 3. The observed brightness profile for L483 at 450 /im 
(upper panel) and 850 /im (lower panel) with the best-fit 
models overplotted (full line). The dashed line indicates 
the beam profile used in the modelling. 



3.3. Results 

In Table ||, the fitted values of the three parameters for 
each source are presented and in Table || the physical pa- 
rameters obtained by scaling according to source distance 
and luminosities are given. As obvious from the \ 2 plots 
in Fig. ||, these parameters have some associated uncer- 
tainties. The value of a is determined typically within 
±0.2 leading to a similar uncertainty in rioo of ±0.2. 
The minimum value of Y gives a corresponding minimum 
value of the outer radius. As discussed above, increasing Y 
only corresponds to adding more material after the outer 
boundary, so if Y is large enough to encompass the point 
where the temperature T reaches 10 K, the radius corre- 
sponding to this temperature can be used as a characteris- 
tic size of the envelope. The region of the envelopes within 
this radius corresponds to the inner 40-50" of the bright- 
ness profiles for all the sources. If one would increase the 
radius further, the typical 120" chop throw for SCUBA 
should be taken into account when comparing the bright- 
ness profiles from the models with the observations. This 
could lead to flatter density profiles with a decreased by 
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Fig. 4. The spectral energy distribution of L483: the 
points indicate the data, the line the best-fit model. 

Table 5. Best fit parameters from DUSTY modelling. 



Source 


y(a) 


a 


T l()0 


L1448-I2 


1800 


1.2 


1.1 


L1448-C 


1600 


1.4 


0.5 


N1333-I2 


900 


1.8 


1.6 


N1333-I4A 


1000 


1.8 


6.5 


N1333-I4B 


1400 


1.3 


0.9 


L1527 


2500 


0.6 


0.1 


VLA1623 


2400 


1.4 


0.7 


L483 


1400 


0.9 


0.3 


L723 


2500 


1.5 


1.0 


L1157 


600 


1.7 


3.4 


CB244 


2600 


1.1 


0.2 


L1489 


1200 


1.8 


0.3 


TMR1 


2000 


1.6 


0.1 


L1551-I5 


1000 


1.8 


1.1 


TMC1A 


1700 


1.9 


1.3 


TMC1 


2900 


1.6 


0.2 


LI 544 


2800 


0.1 


0.1 


L1689B 


3000 


0.1 


0.2 



Notes: ^ Value corresponding to rijri in Table |(| 



0.2 (e.g. Motte fc Andre, 2001). Although the models 



presented here do not extend that far, one should still be 
aware of the possibility that emission is picked up at the 
reference position, which would lead to an overestimate 
of the steepness of the density distribution. The fitted 
brightness profiles and SEDs for all sources are presented 
in Fig. |, § and @. 

The derived power-law indices are for most sources in 
agreement with the predictions from the inside-out col- 
lapse model of a = 1.5 - 2.0. A few YSOs (esp. L1527, 
L483, CB244 and L1448-I2) have density distributions 
that are flatter, but as argued in Sect. 4.2, this can be 
explained by the fact that these sources seem to have sig- 
nificant departures from spherical symmetry. Apart from 
these sources, there is no apparent distinction between the 
class and class I sources in the sample as shown in the 
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Fig. 8. The slope of the density distribution (upper panel) 
and the mass within the 10 K radius (lower panel) vs. 
the bolometric temperature of the sources. Class objects 
are marked by "♦" , class I objects by "A" and pre-stellar 
cores by "■" . VLA1623 and IRAS 16293-2422 have been 
singled out with respectively "•" and "if". The mass of 
the pre-stellar cores in the lower panel is the mass of the 
Bonnor-Ebert sphere adopted for the line modelling. 



plot of power-law slope versus bolometric temperature in 
the upper panel of Fig. |^. There is, however, a signifi- 
cant difference in the masses derived for these types of 
objects, with the class objects having significantly more 
massive envelopes (lower panel of Fig. ||). This is to be 
expected since the bolometric temperature measures the 
redness (or coolness) of the SED, i.e., the amount of en- 
velope material. Finally, there is no dependence on either 
mass or power-law slope with distance, which strengthens 
the validity of the derived parameters. 

3.4. Individual sources 

For some of the sources the derived results are uncertain 
for various reasons, e.g., the interpretation of their sur- 
rounding environment. These cases together with other 
interesting properties of the sources are briefly discussed 
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Source 


Tl 

(AU) 


ri 

(AU) 




rioK 
(AU) 




A^H 2 ,10K 

(cm- 2 ) 


Miok 

(M ) 


n(ri) 
(cm" 




^lOOOAU 

(cm -3 ) 


nioK 
(cm" 




L1448-I2 


7.4 


1 9 
l.O 


X 


10 4 


A £ 

4.0 


X 


10 3 


9 C 
O.O 


X 


i a^3 
1U 


1.5 


O.O 


X 


10 8 


Z.O 


X 


1 A^ 

tu 


A A 

4.U 


X 


10 b 


L1448-C 


9.0 


1.0 


X 


10 4 


0.1 


X 


10 3 


1. / 


X 


1U 


0.93 


£ A 

0.4 


X 


10 8 


1.0 


X 


1 a5 

tu 


A A 

4.U 


X 


10 4 


N1333-I2 


23.4 


z.l 


X 


10 4 


1 o 
l.z 


X 


10 4 


0.0 


X 


1U 


1.7 


1 9 

t.o 


X 


10 9 


t.o 


X 


1 a6 

tu 


1. 1 


X 


10 4 


N1333-I4A 


23.9 


O A 

ZA 


X 


10 4 


A 7 
4. I 


X 


10 3 


o o 
z.z 


X 


to 


2.3 


C A 
O.U 


X 


10 9 


a 9 
6.3 


X 


1 a6 
10 


9 O 

3.8 


X 


10 5 


N1333-I4B 


10.6 


1 c 

t.o 


X 


10 4 


*7 A 
1 .U 


X 


10 3 


9 a 
3.0 


X 


1 a23 

to 


2.0 


D. 1 


X 


10 8 


1 o 
1.0 


X 


1 a6 
10 


1 /I 

1.4 


X 


10 5 


L1527 


4.2 


t.t 


X 


10 4 


O.O 


X 


10 3 


o o 
z.o 


X 


tu 


0.91 


A A 

y.y 


X 


10 6 


9 Q 
O.O 


X 


1 a5 

tu 


l.z 


X 


10 5 


VLA1623 


4.3 


i a 
1.0 


X 


10 4 


9 9 

3.3 


X 


10 3 


O 9 

1.6 


X 


1 n23 

to 


0.22 


t.D 


X 


10 9 


i . ( 


X 


1 A^ 

10 


1.0 


X 


10 5 


L483 


21.5 


3.2 


X 


10 4 


7 o 
1.0 


X 


10 3 


a 9 
9.3 


X 


1 n23 

to 


1.1 


O *7 
Z. ( 


X 


10 9 


1. i 


X 


1 a6 
10 


9 /I 

3.4 


X 


10 4 


L723 


8.2 


z.l 


X 


10 4 


0.4 


X 


10 3 


9 /I 

o.4 


X 


1 a23 

tu 


0.62 


1.4 


X 


10 9 


1.1 


X 


1 a6 

tu 


0.0 


X 


10 4 


L1157 


17.9 


1.1 




10 4 


5.4 




10 3 


1.2 




10 24 


1.6 


3.1 




10 9 


3.3 




10 6 


1.9 




10 5 


CB244 


3.3 


8.7 


X 


10 3 


4.3 


X 


10 3 


6.5 


X 


10 22 


0.28 


2.5 


X 


10 8 


4.8 


X 


10 5 


9.7 


X 


10 4 


L1489 


7.8 


9.4 


X 


10 3 


9.2 


X 


10 3 


1.0 


X 


10 23 


0.097 


7.1 


X 


10 8 


1.2 


X 


10 5 


2.1 


X 


10 3 


TMR1 


6.7 


1.3 


X 


10 4 


1.2 


X 


10 4 


3.5 


X 


10 22 


0.12 


2.1 


X 


10 8 


6.9 


X 


10 4 


1.3 


X 


10 3 


L1551-I5 


24.8 


2.5 


X 


10 4 


1.6 


X 


10 4 


3.8 


X 


10 23 


1.7 


8.2 


X 


10 8 


1.1 


X 


10 6 


7.2 


X 


10 3 


TMC1A 


8.4 


1.4 


X 


10 4 


4.8 


X 


10 3 


4.5 


X 


10 23 


0.13 


3.2 


X 


10 9 


3.7 


X 


10 5 


1.8 


X 


10 4 


TMC1 


3.1 


9.0 


X 


10 3 


4.3 


X 


10 3 


6.9 


X 


10 22 


0.034 


9.0 


X 


10 8 


8.8 


X 


10 4 


8.5 


X 


10 3 


LI 544 


3.0 


8.3 


X 


10 3 


4.0 


X 


10 3 


1.8 


X 


10 22 


0.41 


5.5 


X 


10 5 


3.1 


X 


10 5 


2.7 


X 


10 5 


L1689B 


1.3 


4.0 


X 


10 3 


1.5 


X 


10 3 


2.9 


X 


10 22 


0.096 


2.3 


X 


10 6 


1.2 


X 


10 6 


1.1 


X 


10 6 




Fig. 5. Composite showing the brightness profiles of the sources overplotted with the best fit model at 450 ^im. The 
dashed line indicate the beam profile. 
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Fig. 6. As in Fig. || but for the 850 jtrai data. 



below. The failure of our power-law approach to describe 



the pre-stellar cores will be further discussed in Sect. 4.3 



N1333-I2,-I4: A major factor of uncertainty in the deter- 
mination of the parameters for the sources associated 
with the reflection nebula NGC1333 is the distance 
of these sources. Values ranging from 220 pc (|C erms 
1990) to 350 pc (Herbig fc Jones, 1983) have been sug- 



gested. The latter determination of the distance as- 
sumes that NGC 1333 (and the associated dark cloud 
L1450) is part of the Perseus OB2 association, which 



recent estimates place at 318 ± 27 pc (de Zeeuw et al 
199^). The first value is a more direct estimate based 
on the extinction towards the cloud. N1333-I2 can be 



mo delled uamg 
bee bmea rather large 



cither diat 



ancc — but the 10 K radiua 
i.e., 19000 AU in the caac of 

an assumed distance of 350 pc, while the distance of 
220 pc leads to a radius of 11000 AU that is more 
consistent with those of the other sources. 
Nl 333-14 is one of the best-studied low-mass protostellar 
systems, both with respect to the molecular content 



(Blake et al., 1995) and in interferometric continuum 
studies (Looney et al., 200C). On the largest scales the 



velope, but going to progressively smaller scales shows 
that bo th N1333-I4A and N1333-I4B are multiple in 
nature ( Looney et al., 20"0C| ). The small separation be- 
tween the two sources can cause problems when in- 
terpreting the emission from the envelopes of each of 
the sources. On the other hand the small scale binary 
components of N1333-I4A and N1333-I4B each should 
be embedded in common envelopes and can at most 
introduce a departure from the spherical symmetry. 
L1448-C, -12: The L1448 cloud reveals a complex of 4 or 
more class objects with L1448-C (or L1448-mm) and 
its powerful outflow together wit h the binary protosta r 
L1448-N being well studied (e.g., [Barsony et al., 1998|) 



Also the recently identified L 1448-12 ( O 'Linger et al 



1999) shows typical protostellar properties. The dark 



entire system is seen to be embedded in a single en- 



cloud L1448 itself is a member of the Perseus molec- 
ular cloud complex, for which we adopt a distance of 
220 pc (see above). Note, however, that Ccrnis (1990 ) 
mentions the possibility of a distance gradient across 
the cloud complex so that larger distances may be ap- 
propriate for the L1448 objects. 
VLA 1623: As the first "identified" class object, VLA 
1623 is often discussed as a prototype class object. 
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Fig. 7. Composite showing the SEDs of each source overplotted with the SED from the best fit model. The individual 
SEDs are based on literature searches, with the main ref erences being [Bhirley et al. (2000 ) (class objects and pre- 
stellar cores), ^handler fc Richer (2000| ) (NGC1333-I2), |Sandell et al. (199l|) (NGC1333T4 A,B; 350 and 800 jum), 
IRAS Faint Source Catalog (TMR1, L1489) and 1.3 mm fluxes from [Mottc & Andre (2001| ) for sources included in 
their sample. 



It is, however, not well suited for discussions of the 
properties of these objects because of its location close 



to a number of submillimeter cores (e.g., Wilson et al 
199^). This makes it hard to extract and model the 
properties of this source and might explain why it 
has been claimed t o have a very shall ow density pro- 
file of p rx r" 05 ( Andre et al., 1993) or a const ant 



density outer envelope ( layawardhana et al., 2001 ). If 
the emission in the three quadrants towards the other 
submillimeter cores is blocked out when creating the 
brightness profiles it is found that VLA 1623 can be 
modelled with an almost "standard" density profile 
with a = 1.4, although with rather large uncertain- 
ties. 

L1527: L1527 is remarkable for its rather flat envelope 
profile with a « 0.6. It was one of the sources for 
which phandlcr fc Richer (200C ) estimated the rela- 



tive contribution of the disk and envelope to the total 
flux at 450 and 850 /an and found that the envelope 



contributes by more than 85%, which justifies our use 
of the images at these wavelengths to constrain the en- 



velope. On the other hand, Hogcrheijde et al. (1997) 



estimate the disk contribution at 1.1 mm in a 19" beam 
to be between 30 and 75% of the continuum emission 
(~ 50% for L1527) for a sample of mainly class I ob- 
jects, so it is evident that possible disk emission is a 
factor of uncertainty in the envelope modelling. Disk 
emission would contribute to the fluxes of the inner- 
most points on the brightness profiles and so lead to a 
steeper density profile. 
L723-mm: The most characteristic feature about L723- 
mm is the quadrupolar outflow originating in the cen- 
tral source, which has l ead to the suggestio n that the 
central star is a binary (Girart et al., 1997). 
L483-mm: L483-mm is a good example of a central source 
with an asymmetry yet providing an excellent fit to the 
brightness profile with the simple power-law (see dis- 
cussion in Sect. |4.2|). The source seems to be located in 
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a flattened filament showing up clearly in the SCUBA 
maps (e.g ., |Shirley et al., 2000|) a nd integrated NH 3 
emission flFuller fc Wootten, 2000| ). 
L1157-mm: L1157-mm is not as well known for the pro- 
tostellar source itself as for its bipolar outflow, where 
a large enhancement of chemical species like CH3OH, 
HCN and H2CO is seen (e.g., Bachiller fc Perez| 



Gu ;ierrez, 1997). As a result, the SED of the protostar 
itself is rather poorly determined. Our images in Fig. [| 
show a source slightly extended in the east-west direc- 
tion and with a second object showing up south of the 



source in the 850 /im data. Recently Ghini et al. (2001) 
reported a similar observation of the source and added 
that the southern feature is also seen in the 1.3 mm 
data in the direction of the CO outflow from L1157, 
suggesting an interaction between the outflowing gas 
and the circumstellar dust. 
CB244: CB244 is the only protostar of our sample not 



et 



inc luded in the table of Andre et al. (200C ) . Launhardt 



(1997 ) found that this relatively isolated globule 
indeed has a high submillimeter flux, -Lboi/isubmm > 
2% qualifying it as a class object. It is, however, 
probably close to the boundary between the class 
and I stages: Saraceno et al. (1996] ) found that it falls 
in the area of the class I objects in a Lboi vs. F mm 
diagram. 

L1489: Of the two class I sources in our main sample, 
L1489 has recently drawn attention with the sugges- 
tion that the central star is surrounded by a disk-like 
structure rather than the "usual" envelope for class 



0/1 objects ( Hogerheijde , 
2000). [Hogcrheijde~ ~ 



Sandell, 2000; Hogerheijde 



Sandell examined the SCUBA 
images of this source in comparison with the line emis- 
sion with the purpose of testing the differ ent models 
for the envelope structure - especially the Shu (1977 ) 
infall model. They found that L1489 could not be fitted 
in the inside-out collapse scenario, if both the SCUBA 
images and spectroscopic data are modelled simultane- 
ously. Instead they suggested that L1489 is an object 
undergoing a transition from the class I to II stages, 
revealing a 2000 AU disk, whose velocity structure is 
revealed through high resol ution HCO + interferome- 
ter data ( Hogerheijde, 2001 ). That we actually can fit 
the SCUBA data is neither proving nor disproving this 
result. Hogerheijde fc Sandell in fact remark that it is 
possible to fit the continuum data alone, but that this 
would correspond to an unrealistic high age of this 
source. The modelling of L1489 is slightly complicated 
by a nearby submillimeter condensation - presumably 
a pre-stellar core, which has to be blocked out leading 
to an increase in the uncertainty for the fitting of the 
brightness profile. 
TMR1: This is a more standard class I object showing 
a bipolar nebulosity in the infrared corresponding to 



the outflow cavities of the envelope (Hogerheijde et al. 
199§). 



4. Discussion and comparison 

4.1. Power law or not? 

The first simplistic assumption in the above modelling 
is (as in other recent wor k s, e.g., |5nirlcy et al. (2000D ; 
|Chandler fc Richer (2000| ); |Motte L Andre (200l| ))~that 



the density distribution can be described by a single power 
law. This is not in agreement with even the simplest in- 
fall model, but given the observed brightness profiles of 
the protostars it is tempting to just approximate the den- 
sity distribution wi th a single power law. As an example 
[Shirley et al. (200C ), used this approach citing the results 



of |Adams (1991| ) : if the density distribution can be de- 
scribed by a power-law and the beam can be approximated 
by a gaussian then the outcoming brightness profiles will 
also be a power-law, so a power-law fit to the outer parts of 
the brightness profile will directly reflect the density distri- 
bution. This approach is, however, subject to noise in the 
data and the parts of the brightness profile chosen to be 
considered. On the other hand, our ID modelling clearly 
shows that the data do not warrant more complicated fits 
and that the power-law adequately describes the profiles 
of the sources. Modelling of the detailed line profiles will 
require more sophisticated infall models, since signatures 
for infall exist for around half of the class sources in the 
sample (Andre et al., 2000). However, Schoier et al. (2002[ ) 



show that for the case of IRAS 16293-2422, adopting the 



infall model of Shu (1977) does not improve the quality of 
the fit to the continuum data. 



Chandler & Richer (2000) assumed that the envelopes 
were optically thin. In this case, the temperature pro- 
file can be shown to be a simple power-law as well and 
phandler fc Richer subsequently derived analytical mod- 
els which could be fitted directly to the SCUBA data. For 
the sources included in both samples the derived power- 
law indices agree within the uncertainties. However, as 
illustrated in Fig. ^, the optically thin assumption for the 
temperature distribution is not valid in the inner parts 
of the envelope, especially of the more massive class 
sources, so actual radiative transfer modelling is needed 
to establish the temperature profile, crucial for calcula- 
tions of the molecular excitation and chemical modelling. 

Disk emission can contribute to the fluxes of the in- 
nermost points on the brightness profiles and thus lead 
to a steeper density profile. This is likely to be more im- 
portant in the sources with the less massive envelope, i.e. 
class I objects. In tests where the fluxes within the inner- 
most 15" of the brightness profiles are reduced by 50%, 
the best fit values of a are reduced by 0.1-0.2. This is 
comparable to the uncertainties in the derived value of a, 
but can introduce a systematic error. 

It is interesting to note that there is no clear trend 
in the slope of the density profile with type of object. In 
the framework of the inside-out collapse model, one would 
expect a flattening of the density profiles, approaching 1.5 
as the entire envelope undergoes the collapse. This is not 
seen in the data - actually the average density profile for 
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Fig. 9. The temperature profiles for four selected sources 
with temperature profiles calculated in the Rayleigh- Jeans 
limit for the optically thin assumption overplotted for dif- 
ferent dust opacity laws, k„ oc v& . The dashed line corre- 
sponds to (3 = 1 and the dash-dotted line to (3 — 2. 



the class I objects is slightly steeper than for the class 
objects. On the other hand, the suggestion of an outer 
envelope with a flat density distribution and a significant 



fraction of material as suggested by Jayawardhana ct al 
(200l[j]can also not be confirmed by this modelling. As 
seen from the fits in Fig. [| and ^| the brightness profiles 
only suggest departures from the single power-law fits in 
the outer regions in a few cases, so if such a component 
is present, it is not traced directly by the SCUBA maps. 
The slightly flatter density profiles for the class objects 
could be a manifestation of such an outer component - 
but the density distributions of the sources in this sample 
are typicall y much steeper (i.e. p oc r~ 3 / 2 ) than those 
modelled by Jayawardhana et al. (p oc r -1 / 2 ). 



4.2. Geometrical effects 

As mentioned above, it was assumed in the modelling that 
the sources are spherical symmetric. This may not be en- 
tirely correct considering the structure of YSOs, which 
may be rotating and are permeated by magnetic fields 



leading to polar flattening (Terebey et al., 1984; Galli & 



Shu, 1993) and which ar e definitely associated with molec- 



ular outflow s and jets ( Bachiller fc Tafalla, 1999| ; [Richer 
2000| ). One problem in this discussion is the exact 



et al. 



shape of the error beam of SCUBA. Typically an error 
lobe pickup of 15% in the 850 pm images and 45% in 
the 450 pm maps is estimated. This error beam is not 
completely spherically symmetric, but is also not well- 
established, so in fact ID modelling may be the best that 
can be done using the SCUBA data. 

Myers ct al. (1998) investigated the results of depar- 
ture from spherical symmetry of an envelope when cal- 
culating the bolometric temperature of YSOs seen under 
various inclination angles. With a cavity in the pole re- 



gion, roughly corresponding to the effect of a bipolar out- 
flow, they found that the bolometric temperature could 
increase by a factor 1.3-2.5 for a typical opening angle 
of 25°. A similar line of thought can be applied to our 
modelling: departure from spherical symmetry by having 
a thinner polar region will affect the determination of the 
SED - a source viewed more pole-on will see warmer mate- 
rial which leads to an SED shifted towards shorter wave- 



lengths and accordingly a lower value of tioo- Myers et al 



also argued that the effect on a statistical sample would be 
rather small, e.g., compared to differences in optical depth, 
but for studies of individual sources like in our case, this 
effect might be of importance. 

The brightness profile will also change in the aspherical 
case. In the case of a source viewed edge-on this would re- 
sult in elliptically shaped SCUBA images with the 850 ^m 
data showing a more elongated structure since the smaller 
optical depths material in the polar regions would reveal 
material being warmer and thus having stronger emission 
at 450 pm, compensating for the lack of material in these 
images. If we consider the case where the source is viewed 
entirely pole-on, the image would still appear circular, but 
the brightness profiles would show a steeper increase to- 
wards the center in the 450 pm data (closer to the spher- 
ical case) than the 850 pm data for the same reasons as 
mentioned above. 

The good fits to the brightness profiles given the of- 
ten non-circular nature of the SCUBA images is expected 
based on the mathematical nature of power-law profiles. 
Consider as an example a 2D image of a source described 
by: 



I 



(1) 



where f(9) is a function describing the variation of the 
slope of the brightness profiles extracted in rays along dif- 
ferent directions away from the center position. In t he case 
of a source with a simple density profile p oc r~ p Adams 



(1991 ) showed that such sources will also give images with 



power-law brightness distributions corresponding to the 
case where f(9) = const. As a somewhat simplistic case, 
assume that f(9) is a step function with a power-law slope 
Pi in 9 <E [0, 7r[ and a P2 in 9 e [it, 2ir[. In this case deter- 
mination of the power-law slope will be dominated by the 
flattest of the two slopes: the azimuthal averaged bright- 
ness profile will be 



(I)e oc r~ Pl + r~ 



P2 



(2) 



and due to the power-law decline the term corresponding 
to the shallower slope will dominate the average, especially 
at the larger radii. This has two effects: first the distribu- 
tion of power-law "rays" must be very strongly varying 
in order not to result in a power-law average, and second, 
the density profile for asymmetric sources will be flattened 
compared to more spherical sources. 

As a test, brightness profiles were extracted in angles 
covering respectively the flattest and steepest direction of 
L483 and L1527. Modelling the so-derived brightness pro- 
files gives best-fit density profiles of 0.9 and 1.2, compared 
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to the 0.9 derived as the average over the entire image for 
L483, and 0.6 and 0.8 for L1527 compared to 0.6 from the 
entire image. Thus, the brightness profile and the derived 



e.g., the convolution of the "real" brightness profiles with 
the SCUBA beam. 



For modelling of sources without central heating Evans 



densitj profile might indeed be flattened by the asymme- et al. note that their modelling does not rule out a power- 



try and could account for the somewhat flatter profiles 
found towards some sources. The discrepancies between 
the profiles along different directions are, however, not 
much larger than the uncertainties in the derived power- 
law slope, so these sources could have intrinsic flatter den- 
sity distributions instead; with the present quality of the 
data both interpretations are possible. 



4.3. Pre-stellar cores 

The modelling of the pre-stellar cores is more compli- 
cated than that of the class and I sources, since it is 
not clear whether the cores are undergoing gravitational 
collapse, or are centrally condensed and/or are gravita- 
tionally bound. In the case of thermally supported gravi- 
tationally bound cores, the solution for the density profile 
is the Bo nnor-Ebert sphere flEbert, 1955 ; Bonnor, 1956| ). 
Recently Evans et al. (2001 ) modelled three pre-stellar 
cores (including L1689B and L1544 in our sample) and 
found t hat they coul d be well fitted by Bonnor-Ebert 
spheres. |Evans et al. also found that the denser cores, 
L1689B and L1544, were those showing spectroscopic signs 
of contraction thus suggesting an evolutionary sequence 
with LI 544 as the pre-stellar cores closest to the collapse 
phase. This is supported as well by millimeter observa- 
tions of this core which show a dense inner region ( Tafalla 
ct al., 1998; ; |Ward- Thompson et al., 1999| ). 



Evans et al. also discussed the properties of non- 
isothermal models and found that the quality of the fit to 



the data is similar to the isothermal case. Ward-Thompson 



et al. (2002) examined ISOPHOT 200 ^m data for a sam- 



ple of pre-stellar cores (including L1689B and L1544) and 
found that none of the cores had a central peak in tem- 
perature and that they could all be interpreted as be- 
ing isothermal or having a temperature gradient with a 
cold center as res ult of external heatin g by the interstel- 
lar radiation field, ^ucconi ct al. (2001 ) derived analytical 
formulae for the dust temperature distributions in pre- 
stellar cores showing that these cores should have temper- 
atures varying from 8 K in the center to around 15 K at 
the boundary. These equations will be useful for more de- 
tailed modelling of the continuum and line data, but for 
the present purpose the isothermal models are sufficient. 

Modelling the pre-stellar cores using our method is not 
possible as is illustrated by the best fits for the pre-stellar 
cores shown in Fig. ^|-|7]. Given the observational evidence 
that these cores do not have central source of heating, 
what is implicitly assumed in the DUSTY modelling, it 
is on the other hand comforting that our method indeed 
distinguishes between these pre-stellar cores and the class 
0/1 sources and that the obtained fits to the brightness 
profiles of the latter sources are not just the results of, 



law envelope density distribution for L1544, as opposed to, 
e.g., the case for L1689B: this would imply an evolution- 
ary trend of the pre-stellar cores having a Bonnor-Ebert 
density distribution, which would then evolve towards a 
power-law density distribution with an increasing slope 
as the collapse progresses. In the modelling of spectral 
line data for the pre-stellar cores we adopt an isothermal 
Bonnor-Ebert sphere with n c ~ 10 6 for both L1689B and 
L1544 as this was the best fitted isothermal model in the 
work of Evans et al. (200l| ). 



5. Monte Carlo modelling of CO lines 

5.1. Method 

One main goal of our work is to use the derived physical 
models as input for modelling the chemical abundances 
of the various molecules in the envelopes. To demonstrate 
this approach, modelling of the first few molecules, C 18 
and C 17 0, is presented here. This modelling also serves 
as a test of the trustworthiness of the physical models: is 
it possible to reproduce realistic abundances for the mod- 
elled molecules? 



The ID Monte Carlo code developed by |Hogcrheijde fe 
van dcr Tak (200C| ) is used t ogether with the revised colli- 
sional rate coefficients from Flower (2001 ) , and a constant 
fractional abundance over the entire range of the envelope 
is assumed as a first approximation. Furthermore the dust 
and gas temperatures are assumed identical over the en- 
tire envelope and any systematic velocity field is neglected. 
In the outer parts of the envelopes, the coupling between 



gas and dust may break down (e.g. Ceccarclli ct al., 1996 



Doty fc Neufeld, 1997) leading to differences between gas 



and dust temperatures of up to a factor of 2. 

With the given physical model and assumed molecular 
properties, there are two free parameters which can be ad- 
justed by minimizing x 2 to model the line profiles for each 
molecular transition: the fractional abundance [X/H2] and 
the turbulent line width Vd- Since emission from the am- 
bient molecular cloud might contribute to the lower lying 
1-0 transitions, the derived abundances are based only on 
fits to the 2-1 and 3-2 lines. The observed and modelled 
line intensities are summarized in Tables 0-^, whereas the 
parameters for the best fit models are given in Table [n| 
While a turbulent line width of 0.5 — 1.0 km s _1 is needed 
for most sources to fit the actual line profile, the modelled 
line strengths are only weakly dependent on this param- 
eter compared to the fractional abundance. For example, 
for L723 one derives C ls O abundances between 3.8 x 10 -8 
and 4.0 x 10~ 8 from respectively V D = 1.1 km s" 1 (FWHM 
~ 1.9 km s- 1 ) to Vb = 0.7 (FWHM ~ 1.2 km s" 1 ), 
which should be compared to the [C 18 0/H 2 ]= 3.9 x 10~ 8 
and Vb = 0.8 km s" 1 (FWHM ~ 1.4 km s" 1 ) quoted 
in Table |ll|. The C 18 and C 17 spectra for L723 and 
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Table 8. As in Table but for the IRAM 30m observa- 
tions - all from November 2001 observing run. 



V km s 

LSR 



V | km s" 

LSR 



5 
4 

g 3 

I 2 

1 


2.5 
2.0 

£ 1 - 5 
| 1 .0 

i — 

. 5 
0.0 



i 


C ,8 (2-1) 




C ,8 (3-2) j 




fc ,7 (2-1) - 




C 17 (3-2) i 

\ 



5 10 15 5 10 15 



V [km s 

LSR L 



V [km s 

LSR L 



Fig. 10. The C 18 and C 17 spectra for L723 (upper four 
panels) and N1333-I2 (lower four panels) overplotted with 
the best fit models from the Monte Carlo modelling. 

N1333-I2 with the best fit models overplotted are shown 
in Fig. [To[ The revised collisional rate coefficients from 



Flower (2001 ) adopted for this modelling typically change 
the derived abundances by less than 5% compared to re- 
sults obtained from simulations with previously published 
molecular data. The relative intensities of the various tran- 
sitions and so the quality of the fits remain the same. 

As it is evident from Tables ffiH^ the 2-1 and 3-2 lines 
can be well modelled using the above approach, while the 

1- lines are significantly underestimated from the mod- 
elling, especially in the larger Onsala 20m beam. This in- 
dicates that the molecular cloud material may contribute 
significantly to the observations or that the assumed outer 
radius is too small. The importance of the latter effect was 
tested by increasing the outer radius of up to a factor 2.5 
for a few sources and it was found that while the 1-0 and 

2- 1 line intensities could increase in some cases by up to a 
factor of 2, the 3-2 line intensities varied by less than 10%, 
illustrating that the 3-2 line mainly trace the warmer (> 
30 K) envelope material. 

The importance of the size of the inner radius was 
tested as well. For L723 fixing the inner radius at 50 AU 
rather than 8 AU increases the best fit abundance by ~ 



Source 


C iv 


O 1-0 


C iv 


O 2-1 




Obs 


Mod 


Obs 


Mod 


L1448-C 


0.58 


0.71 


2.2 


2.3 


N 1333-12 


1.1 


0.61 


3.4 


2.5 


N1333-I4A 


0.99 


0.54 


2.9 


2.3 


N1333-I4B 


0.73 


0.47 


1.8 


1.5 


L1527 


1.1 


0.90 


2.4 


2.0 


L1489 


0.37 


0.27 


1.1 


1.1 


TMR1 


0.69 


0.37 


1.4 


1.3 



Notes: The intensities are the total observed intensity summed 
over the hyperfine-structure lines. 

Table 9. As in Table @ and § for the Onsala 20m obser- 
vations - all from March 2002 observing run. 



Source 


C 18 


O 1-0 


c rv 


O 1-0 




Obs 


Mod 


Obs 


Mod 


L1448-I2 


3.2 


0.5 


1.1 


0.2 


L1448-C 


3.6 


1.3 




0.5 


N1333-I2 


5.6 


1.4 


1.6 


0.5 


N1333-I4A 


4.2 


1.9 




0.4 


N1333-I4B 


4.1 


1.5 




0.4 


L1527 


2.4 


1.5 




0.8 


L483 


3.8 


1.7 


2.2 


0.6 


L723 


1.6 


0.6 


0.3 


0.1 


L1157 


1.3 


0.3 


0.2 


0.1 


CB244 


1.9 


1.1 




0.2 


L1489 


1.7 


1.0 




0.2 


TMR1 


1.7 


1.5 




0.3 


LI 544 


1.7 


0.3 


0.6 


0.2 



Notes: '"'The intensities are the total observed intensity 
summed over the hyperfine-structure lines. 



5% to [C 18 0/H 2 ]= 4.1 x 10~ 8 without changing the qual- 
ity of the fit significantly (drop in \ 2 01 ~ 0.1). This is 
also found for other sources in the sample and simply il- 
lustrates that increasing the inner radius corresponds to a 
(small) decrease in mass of the envelope, so that a higher 
CO abundance is required to give the same CO intensities. 
Yet, it is good to keep in mind that the inner radius of 
the envelope may be different, and even though its value 
does not change the results for CO, it might for other 
molecules, e.g., CH3OH and H2CO, which trace the inner- 
warm and dense region of the envelope. 

5.2. CO abundances 

The fitted abundances are summarized in Table [H] and 
plotted against the envelope mass of each individual 
source m Fig. [tl]. As can be seen from the values of the 
reduced \ 2 f° r the fit to the data for each isotope, the 
model reproduces the excitation of the individual species. 
Even in the "worst case" of N1333-I2, the lines do indeed 
seem to be well- fitted by the model (see Fig. [l(]). One 
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Table 7. The observed C 18 and C 17 lines from the JCMT observations compared to the modelled intensities. 



Source 








(2-1) 


(3-2) 




(kms" 1 )*"' 


(km s 


- 1 ) Obs (b) Mod (c) 


Obs (b) 


Mod (c) 


C 18 


L1448-I2 


0.7 


4.5 




1.6 


1.9* 


1.9 


L1448-C 


1.2 


5.3 


4.3 


3.3 


3.3* 


3.9 


N1333-I2 


1.4 


7.7 


5.8 


3.9 


4.7* 


5.7 


N1333-I4A 


1.7 


7.2 


5.5 


5.4 




7.1 


N1333-I4B 


2.1 


7.4 


4.3 


3.8 


4.0* 


4.4 


L1527 


0.7 


5.9 


4.4 


2.9 


2.3 


2.7 


VLA1623 


1.0 


3.5 


12.1: 


12.1 


12.6: 


12.6 


L483 


1.0 


5.2 


4.1 


3.7 


3.6 


3.8 


L723 


1.6 


11.2 


2.2 


1.9 


2.2* 


2.4 


L1157 


0.8 


2.6 


1.7 


1.2 


1.4 


1.7 


CB244 


1.1 


4.4 


3.3 


2.9 


2.9 


3.1 


L1489 


2.1 


7.2 


2.7 


2.6 


3.8 


3.9 


TMR1 


1.5 


6.3 


4.0 


3.3 


3.8 


4.3 


L1551-I5 


2.1 


7.2 


7.1 


6.1 


8.7 


9.7 


TMC1A 


1.5 


6.6 


1.3 


1.5 


2.8 


2.3 


TMC1 


1.5 


5.2 


2.3 


2.7 


4.3 


3.2 


L1544 


0.3 


7.6 




0.55 


0.30 


0.30 


L1689B 


0.6 


3.6 


3.6 


3.3 


1.8 


1.9 


C 17 


L1448-I2 


0.9 


4.0 




0.61 


0.71* 


0.72 


L1448-C 


1.4 


5.0 


1.5 


1.4 


1.7* 


1.7 


N1333-I2 


1.3 


7.5 




1.4 


1.8* 


2.1 


N1333-I4A 


1.3 


6.7 




1.1 


1.4* 


1.6 


N1333-I4B 


1.4 


6.8 




0.9 




1.1 


L1527 


0.5 


6.0 


1.9 


1.7 


1.3 


1.6 


VLA1623 


0.8 


3.7 




3.8 


4.6:* 


4.5 


L483 


0.8 


5.3 


1.5 


1.3 


1.2* 


1.3 


L723 


1.3 


11.0 


0.59 


0.43 


0.49* 


0.57 


L1157 


1.0 


2.7 


0.51 


0.47 


0.65* 


0.69 


CB244 


1.0 


4.2 


0.94 


0.68 


0.60* 


0.73 


L1489 


3.0 


7.3 




0.56 


0.86* 


0.88 


TMR1 


1.6 


6.0 




0.69 


0.84* 


0.94 


L1551-I5 


2.1 


7.2 


2.2 


1.8 


2.6 


3.1 


TMC1A 


1.5 


6.6 




0.45 




0.73 


TMC1 


1.5 


5.2 




0.83 




1.0 


L1544 


0.6 


7.5 


0.26 


0.26 


*« 


0.14 


L1689B 


0.6 


3.7 


0.71 


0.85 


0.59* 


0.44 



Notes: (a) The width (FWHM) of the lines. (6) Observed intensities, / T mb dV in K km s _1 : our own observations are marked with 
'*' and lines where double gaussians were fitted and the broadest component subtracted marked with ':'. '^Modelled intensities 
in K km s _1 . ^C ir O 3-2 not detected towards L1544 in integrations corresponding to an RMS of 0.1 K (T£)- 



source, VLA1623 shows a remarkably high ratio between 
the C ls O and C 17 abundances of 12.4 and high C 18 
and C 17 abundances in comparison to the rest of the 
class sources. Given the location of VLA 1623 close to 
a dense ridge of material and the associated uncertainties 
in the physical models of this source, it may reflect prob- 



lems iii the inudel laLhei Llia.ii being a leal proper l 
the source. Fui tiro lesL uf Lliu suuices, Lhu laLiu between 

the C 18 and C 17 abundances is found to be 3.9 ± 1.3 
in agreement with the expected value from, e.g., the lo- 



Lly uf 



cal interstellar medium of 3.65 (Pcnzias, 1981; Wilson & 



Rood, 1994) - another sign that the model reproduces the 
physical structure of the envelopes and that no systematic 



calibration errors are introduced by using data from the 
various telescopes and receivers. 

It is evident that the class objects (except VLA 
1623) and pre-stellar cores show a high degree of deple- 
tion compared to the expected abundances of [C 18 0/H2] 
of 1.7 x 10~ 7 from nearby dark clouds (Frerking et al. 
79821 ) and [C 17 0/H 2 ] of 4.7 x 1CT 8 assuming 18 0: 17 of 
3.65. With this isotope ratio and assuming 16 0: 18 equal 



to 540 (Wilson & Rood, 1994), the average abundance 
for the class I sources is (1.1 ± 0.9) x 10 -4 and for the 
class sources and pre-stellar cores (2.0 ± 1.3) x 10~ 5 . 
The error bars illustrate the source to source variations 
and uncertainties in classifying borderline class 0/1 ob- 
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Table 10. As Tables for CSO measurements from the 
literature. 



Source 


C ls O 3-2 (a) 


C 17 3-2 (a) 




Obs Mod 


Obs Mod 


N 1333-14 A 


4.9 4.9 


1.0 1.0 


N1333-I4B 


3.3 3.3 


0.80 0.80 


C 17 2-l< 6 > 


TMCIA 


0.28 0.27 




TMCl 


0.53 0.51 





Notes: (a) Fr om |Blakc et al. (199^ ), half-power beam width of 
20". (b) From|Ladd et al. (1998|), half-power beam width of 33". 



jects li ke CB 244, L1527 and L1551-I5. Previously pasclli 
et al. (1999) derived the C 17 abundance for one of the 



pre-stellar cores, L1544, and our abundance agrees with 
their estimate within the uncertainties. It is interesting to 
see that the class I objects have higher CO abundances 
close to the molecular cloud values, indicating that the 
class objects indeed seem to be closer relate d to the pre- 
stellar cores in this sense. Van der Tak et al. ( 2000b ) like- 
wise found a trend of increasing CO abundance with mass- 
weighted temperature for a sample of high-mass YSOs and 
suggested that this trend was due to freeze-out of CO in 
the cold objects. 

The derived abundances are uncertain due to several 
factors. First, the physical model and its simplicity and 
flaws as discussed above. Second, the assumed dust opac- 
ities will affect the results: the dust opacities may change 
with varying density and temperature in the envelope, and 
will depend on the amount of coagulation and formation 
of ice mantles. These effects tend to increase k v with in- 
creasing densities or lower temperatures, which will lower 
the derived mass and thus increase the abundances nec- 
essary to reproduce the same line intensities. Comparing 
the models of dust opacities in environments of different 



densities and types of ice mantles as given by Osscnkopf 



& Hen: ling (1994 ) indicates, however, that such variations 
should be less than a factor of two and so cannot explain 
the differences in the derived CO abundances. 

The unknown contribution to the lower lines from e.g. 
the ambient molecular cloud may affect the interpretation. 
Even for the 2-1 transitions the cloud material may con- 
tribute, leading to higher line intensities than predicted 
from the modelling. Indeed as judged from Tables |^-pi3| 
there is a trend that while the intensities of 2-1 lines are 
underestimated by the modelling, the 3-2 lines are overes- 
timated. This then means that the derived CO abundances 
are upper limits - at least for the warmer regions traced 
by the 3-2 lines. 

As noted above the gas temperature could be lower by 
up to a factor of two in the outer parts of the envelope 



( |Ceccarelli et al., 1996| ; poty fc Ncufcld, 1997| ), but this 
again would mainly affect the lines tracing the outer cold 
regions i.e. the 1-0 and 2-1 lines compared to the 3-2 
lines. At the same time the gas temperatures from these 
detailed models may not be appropriate for our sample. 



Our sources all have lower luminosities than modelled e.g. 
by Ceccarelli et al. (1996), leading to envelopes that are 
colder on average, so that the relative effects of the de- 
coupling of the gas and dust are smaller. Another effect 
is external heating, which as discussed in Sect. 4.3 may 
lead to an increase in the dust and gas temperatures in 
the outer parts. 

A point of concern is also whether steeper density 
gradients could change the inferred abundances in light 
of the discussion about geometrical effects (Sect. 4.2). 
Steepening of the density distribution would tend to shift 
material closer to the center and so towards higher temper- 
atures. This would correspondingly change the ratio be- 
tween the 2-1 and 3-2 lines towards lower values (stronger 
3-2 lines) and more so for the C 18 data than the C 17 
data, since it traces the outer less dense parts of the enve- 
lope. Altogether none of the effects considered change the 
conclusion that the abundances in the class objects are 
lower than those found for the class I objects. 

One important conclusion of the derived abundances 
is that one should be careful when using CO isotopes to 
derive the H2 mass of, e.g., envelopes around young stars 
assuming a standard abundance. With depletion the de- 
rived H2 envelope masses will be underestimated. Another 
often encountered assumption, which may introduce sys- 
tematic errors, is that the lower levels are thermalised and 
that the Boltzmann distribution can be used to calculate 
the excitation and thus column density of a given molec- 
ular species. In Fig. [l^, the level populations for C 18 
in the outer shell of the model of two sources, L723 and 
TMR1, are shown together with the variation of the ratios 
of the level populations from the Monte Carlo modelling 
by assuming LTE. For the more dense envelope around the 
typical class object, L723, the LTE assumption gives ac- 
curate results within 5-10% in the lower levels, but some- 
what more uncertain results for the higher levels. For the 
less dense envelope around TMR1 (class I object) it is, 
however, clear that the levels are subthcrmally excited and 
the LTE approximation provides a poor representation of 
the envelope structure. 

5.3. CO abundance jump or not? 

The large difference in CO abundance found between the 
class sources and pre-stellar cores on the one side and 
the class I objects on the other side warrants further dis- 
cussion. 

The apparently low CO abundances and the possible 
relation to freeze-out of CO raises the question whether 
the assumption of a constant fractional abundance is real- 
istic: freezing out of pure CO-ice and isotopes is expected 
to occur at roughly 20 K under interstellar conditions 
(e.g., [Bandford fc Allamandola, 1993 ), so one would expect 
to find a drastic drop in CO abundance in the outer parts 
of the envelope. Due to the uncertainties in the properties 
of the exterior regions, however, a change in abundance at 
20 K can neither be confirmed nor ruled out. As Table [?] 
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Table 11. CO abundances using the best models of Table ^. 



Source 




[C ls O/H 2 ] 


x 2 d 


[C 17 0/H 2 ] 


i 


r« /"0 


[CO/H 2 ] lc) 


L1448-I2 


0.5 


l.OxlO -8 




3.5xl0 -9 




2.9 


6.1xl0~ b 


L1448-C 


0.7 


5.6xl0 -8 


2.1 


2.2xl0~ 8 


0.07 


2.5 


3.7xl0 -5 


N1333-I2 


0.8 


4.1xl0" 8 


3.9 


1.3xl0~ 8 


2.5 


3.2 


2.4xl0~ 5 


N1333-I4A 


0.7 


1.9xl0 -8 


< 0.01 


2.8xl0~ 9 


0.7 


6.8 


7.9xl0~ 6 


N1333-I4B 


0.6 


2.8xl0" 8 


0.30 


5.9xl0~ 9 


< 0.01 


4.7 


1.3 xlO -5 


L1527 


0.4 


4.9xl0" 8 


3.5 


2.6xl0" 8 


1.0 


1.8 


3.9xl0" 5 


VLA1623 


0.6 


1.0xl0~ 6 


< 0.01 


8.1xl0~ 8 




12.3 


1.6xl0~ 4 (e) 


L483 


0.6 


2.5xl0" 8 


0.28 


7.8xl0~ 9 


1.1 


3.2 


1.4xl0 -5 


L723 


0.8 


3.9xl0" 8 


0.83 


8.1xl0~ 9 


2.5 


4.8 


1.9 xl0~ 5 


L1157 


0.6 


l.OxlO" 8 


3.3 


3.6xl0~ 9 


0.23 


2.8 


6.2xl0~ 6 


CB244 


0.6 


8.0xl0~ 8 


0.47 


1.6xl0~ 8 


3.1 


5.0 


3.7xl0~ 5 


L1489 


1.2 


2.2xl0~ 7 


0.13 


4.5xl0~ 8 


0.04 


4.6 


l.OxlO" 4 


TMR1 


0.9 


4.5xl0~ 7 


1.4 


8.6xl0~ 8 


1.0 


4.5 


2.0xl0~ 4 


L1551-15 


0.9 


5.6x 10 -8 


0.80 


1.5x 10 -8 


1.3 


3.7 


O A w 1 A — 5 

3. Ox 10 


TMC1A 


0.7 


4.3xl0" 8 


1.3 


1.2xl0~ 8 




3.6 


2.3 xl0~ 5 


TMC1 


0.7 


3.6xl0" 7 


2.1 


1.0xl0~ 7 




3.6 


2.0xl0" 4 


L1544 


0.3 


6.8xl0~ 9 




3.1xl0~ 9 




2.2 


4.9xl0~ 6 


L1689B 


0.5 


5.1xl0" 8 


0.36 


1.0xl0~ 8 


2.45 


5.1 


2.4xl0~ 5 


IRAS16293-2422 (/) 




6.2xl0~ 8 




1.6xl0~ 8 




3.9 


3.3xl0~ 5 



Notes: 



indicate abundances where only one line where available to constrain the fit. '"^ Turbulent velocity in km 



The 18 0/ 17 isotope ratio (or [C 18 0/H 2 ]/[C 17 0/H 2 ]) (c) Derived CO main isotope abundance averag ed from the C ls O and 
C 17 measurements assuming ls O/ 17 of 3.65 and 16 0/ ls O of 540 ( [Penzias, 198l[ |Wilson fc Rood, 1994] ) (e) Based on the C 17 
measurements only. (/) IRAS 16293-2422 included for comparison; for details see [Schoier et al. (2002|). 




l°g(n(H,)) 



Fig. 12. The derived level populations for C 18 in the 
envelopes around the class object L723 (left column) and 
class I object TMR1 (right column). In the upper panels 
the level populations from the modelling in the outer shell 
of the envelope are shown (dashed line) together with the 
predictions from the LTE assumptions (solid line) . In the 
lower panels, the ratio of the resulting level populations 
from the modelling and the LTE predictions are shown 
with varying density for the J = to 3 levels. 



indicates both the intensities of the 2-1 and 3-2 lines of 
C 18 and C 17 can be fitted with a constant fractional 
abundance for most sources. One can introduce strong de- 
pletion through a "jump" in the fractional abundance in 



the region of the envelope with temperatures lower than 
20 K. This naturally leads to lower line intensities of espe- 
cially the 2-1 and 1-0 lines, but also the 3-2 lines. Since 
it is mainly the 3-2 line that constrains the abundance in 
the inner part, it is possible to raise the abundances of the 
warm regions by up to a factor 2, if CO in the outer part 
is depleted by a factor of 10 or more. The modelled 1-0 
and 2-1 line intensities then, however, also become weaker, 
which one has to compensate for by introducing even more 
cold (depleted) material outside the 10 K boundary, ac- 
counting for 50% or more of the observed 2-1 line emission 
and almost all the 1-0 emission. 

If the dust opacity law varies with radial distance in 
the envelope, increasing k„ with the higher densities, the 
CO abundances in the warmer regions could be higher 
by a factor of two. Combining these effects would then 
lower the derived CO evaporation temperature in the en- 
velope material towards the 20 K evaporation temperature 
of pure CO ice. A changing dust opacity law would be an 
interesting result in itself, which possibly can be confirmed 
or disproved through modelling of other chemical species 
and/or by comparing the exact line profiles with realistic 
models of the velocity field. Without such an effect in- 
troduced, however, the results presented here favors the 
somewhat higher evaporation temperature for the CO. 

The differences between the class and I objects and a 
warmer evaporation temperature may be consistent with 
new laboratory experiments on the trapping and evapo- 
ration of CO on H 2 ice by pollings et al. (2002|). Their 



experiments refer to amorphous H2O ice accreted layer- 
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Fig. 11. The fitted abundances vs. envelope mass (left) and bolometric temperature (right) of each source for re- 
spectively the C 18 data (top) and C 17 data (bottom). The sources have been split into groups (class 0, class I 
and pre-stellar cores) with VLA 1623 and IRAS 16293-2422 separated out using the same symbols as in Fig. ||: class 
objects are marked with class I objects with "A", pre-stellar cores with VLA 1623 with "•" a nd IRAS 
16293- 2422 with "-^r" . The vertical lines in the figures illustrate the abundances in quiescent dark clouds from Frerking 
et al. ([L982D. 



by-layer so that it has a porous ice structure. This situa- 
tion may be representative of the growth of H2O ice layers 
in pre-stellar cores and YSO envelopes. CO is deposited 
on top of the H2O ice. When the sample is heated from 
10 to ~30 K (laboratory temperatures), some of the CO 
evaporates, but another fraction diffuses into the H2O ice 
pores. Heating to 30-70 K allows some CO to desorb, but 
a restructuring of the H2O ice seals off the pores, and the 
remaining CO stays trapped until at least 140 K. Under in- 
terstellar conditions, the temperatures for these processes 
may be somewhat lower, but it does indicate that a sig- 
nificant fraction of CO can be trapped in a porous surface 
and that evaporation may occur more gradually from 30 K 
to ^60 K. A similar property was suggested by Ceccarclli 
et al. ( fe00l| ), who found that the dust mantles in the enve- 
lope around IRAS 16293-2422 had an onion-like structure 
with H2CO being trapped in CO rich ices in the outermost 
regions and with the ices becoming increasingly more H2O 
rich when moving inwards toward higher temperatures. In 



this scenario, it is not surprising that even the 3-2 lines 
tracing the warmer material indicate low CO abundances. 
Observations of even higher lying CO rotational lines (e.g., 
6-5) would be needed probe the full extent of this evapo- 
ration. Also infrared spectroscopy of CO ices may reveal 
differences between the class and I objects. 



6. Conclusion 

This is the first paper in a survey of the physical and chem- 
ical properties of a sample of low-mass protostellar ob- 
jects. The continuum emission from the envelopes around 
these objects has been modelled using the ID radiative 
transfer code, DUSTY, solving for the temperature distri- 
bution assuming simple power-law density distributions of 
the type p oc r~ Q . For the class and I objects in the sam- 
ple, the brightness profiles from SCUBA 450 and 850 /im 
data and the SEDs from various literature studies can be 
successfully modelled using this approach with a in the 
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range from 1.3 to 1.9 within ^10000 AU with a typical 
uncertainty of ±0.2, while it fails for the pre-stellar cores. 
For four sources the profiles are indeed flatter than pre- 
dicted by, e.g., the models of [Shu (1977 ), but it is argued 
that this could be due to source asymmetries and/or the 
presence of extended cloud material. Taking this into ac- 
count, no significant difference seems to exist between the 
class and I sources in the sample with respect to the 
shape of the density distribution, while, as expected, the 
class objects are surrounded by more massive envelopes. 

The physical models derived using this method have 
been applied in Monte Carlo modelling of C 18 and C 17 
data, adopting an isothermal Bonnor-Ebert sphere as a 
physical model for the pre-stellar cores. The 2-1 and 3-2 
lines can be modelled for all sources with constant frac- 
tional abundances of the isotopes with respect to H 2 and 
an isotope ratio [ 18 0/ 17 0] of 3.9, in agreement with the 
"standard" value for the local interstellar medium. The 
1-0 lines intensities, however, are significantly underesti- 
mated in the models compared to the observations, in- 
dicating that ambient cloud emission contributes signifi- 
cantly or that the outer parts of the envelopes are not well 
accounted for by the models. The derived abundances in- 
crease with decreasing envelope mass - with an average 
CO abundance of 2.0 x 10~ 5 for the class objects and 
pre-stellar cores, and 1.1 x 10~ 4 for the class I objects. 
The 3-2 lines indicate that the lower CO abundance in 
class objects also applies to the regions of the envelopes 
with temperatures higher than ^20 — 25 K, the freeze- 
out of pure CO ice. This feature can be explained if a 
significant fraction of the solid CO is bound in a (porous) 
ice mixture from where it does not readily evaporate. The 
physical models presented here will form the basis for fur- 
ther chemical modelling of these sources. 
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